Intermittence and roughening of periodic elastic media 
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We analyze intermittence and roughening of an elastic interface or domain wall pinned in a 
periodic potential, in the presence of random-bond disorder in (1+1) and (2+1) dimensions. Though 
the ensemble average behavior is smooth, the typical behavior of a large sample is intermittent, and 
does not self-average to a smooth behavior. Instead, large fluctuations occur in the mean location 
of the interface and the onset of interface roughening is via an extensive fluctuation which leads 
to a jump in the roughness of order A, the period of the potential. Analytical arguments based on 
extreme statistics are given for the number of the minima of the periodicity visited by the interface 
and for the roughening cross-over, which is confirmed by extensive exact ground state calculations. 
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I. INTRODUCTION 



The properties of extended, elastic manifolds like do- 
main walls in magnets or contact lines of liquids on solid 
substrates becomes very varied if one introduces some 
disorder. Defects on a surface or impurities in a mag- 
net often pin such interfaces. The recent interest in their 
physics follows from the observation that the energetics 
in the presence of randomness is obtained by optimiz- 
ing the configuration of the manifold Q . A competition 
between elasticity and the random potential arises. It re- 
sults in scale-invariance described by the roughness expo- 
nent that measures the geometrical fluctuations, and the 
energy fluctuation exponent that measures the variation 
of the manifold energy around its mean. It is also related 
to the energy scales of excitations from the state of mini- 
mum energy. The experimental interest in these systems 
arises in particular due to the energetics: time-dependent 
phenomena like creep and coarsening (in magnets) follow 
slow, activated dynamics dictated by the energy barriers 
that can be described with such exponents 0. 

Frequently manifolds also experience a periodic poten- 
tial. In the case of superconductors, one periodicity is 
due to the rotational invariance of the phase. A sec- 
ond periodicity is induced when flux lines form a lattice. 
Similarly, in the case of charge density waves (CDW) or 
domain walls in magnets, one periodicity is due to the 
underlying lattice structure, and a second is due to the 
self-organized periodicity of the CDW or magnetic do- 
mains themselves. Generic models for these phenomena 
are called periodic elastic media (PEM) and are the focus 
of this work. As noted recently the asymptotic behavior 
of the PEM class depends on the type of periodicity, with 
the case of a periodic surface tension being in one uni- 
versality class ||, while the case of an applied periodic 
potential is in another [Q . In this work we are interested 
in the case of an applied periodic potential, in particular 
the intermittent behavior of interfaces which experience 
a competition between pinning due to the periodic po- 



tential and pinning due to random bond disorder. 

The paper is organized as follows. Sec. || introduces 
the Hamiltonian of periodic elastic media and describes 
two intermittent behaviors involved in PEM, when the 
amplitude of the applied periodicity is changed. The Sec- 
tion also includes a discussion of the numerical method 
used. In Sec. Ill the first type of the intermittent be- 
havior, jumps of manifolds, is analyzed using extremal 
statistics and is demonstrated with numerical simula- 
tions. Sec. IV discusses the second type of the inter- 



mittent behavior, the roughening of the manifolds, with 
the aid of droplet arguments and further numerics. In 
Sec. [v] the roughening behavior is studied in {10} and 
{100} oriented lattices which have a lattice induced peri- 
odicity; we compare these systems with other PEM. The 



paper is finished with conclusions in Sec. VI 



II. PERIODIC ELASTIC MEDIA 

The continuum Hamiltonian that describes the com- 
petition between elasticity, periodicity and randomness 
is given by, 
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{Vh(r}} 2 + r 1 {h(r)} + V p {h(r)} dr. (1) 



Here h(r) is a single valued height variable, and r is a 
[d — l)-dimensional vector. V p is a periodic potential (of 
amplitude Vq and wavelength A) in the height direction 
and r/{h(r)} is the disorder, which we take to be of the 
random bond type with delta-function correlations. The 
physics of manifolds described by Eq. (Q) has been dis- 
cussed recently since there may exist a roughening tran- 
sition that separates an algebraically rough regime from 
a logarithmically rough one as the potential strength is 
varied Q]. However, in the dimensions considered here 
[d = (1 + 1), (2 + 1)] these manifolds are always rough at 
large enough length scales ||. The issues we raise here 
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arise in all dimensions, and so we numerically illustrate 
them in (1 + 1)- and (2 + l)-dimensional systems. 

We calculate the exact location and morphology of in- 
terface ground-states for a given configuration of bond 
disorder. For this configuration of bond disorder we 
vary the amplitude of the periodic modulation. Interfaces 
which experience this combination of a periodic potential 
and random bond disorder show a variety of intermittent 
behavior as the amplitude of the potential, V , is varied. 
Two types of intermittence which we study in detail are: 
intermittent jumps in the center of mass location of the 
interface; and intermittent jumps in the roughness of the 
interface. The first type is most easily discussed at strong 
pinning (large values of the key ratio v = Vq\J / 5 J) where 
the interface is always pinned near a minimum of the pe- 
riodic potential, but it jumps between different minima 
as v is varied. It does this to maximize the energy gain 
due to small fluctuations about a flat interface. In the 
limit of large system sizes there can be an infinite num- 
ber of such jumps with, of course, no overlap between the 
ground states of interfaces in different minima. We de- 
velop a scaling theory to demonstrate that the number of 
minima explored as v is varied over a finite range is of or- 
der h\(Lh)i where Lh is the system size in the /i-direction 
in which the manifold fluctuates. Such intermittence is 
similar to the chaos seen in spin glasses (where it implies 
vanishing overlap between spin configurations) and is re- 
lated to replica symmetry breaking || . It is also a close 
cousin of the phenomenon that takes place if the disorder 
is changed randomly j?J . 

A second type of intermittence occurs when it becomes 
energetically favorable to form a large domain excitation. 
This means that a finite fraction of the interface is in one 
minimum of the potential, while another finite fraction 
is in an adjacent minimum. These large fluctuations are 
the classical "Imry-Ma" -type droplets and have a linear 
extension of the order of the system size. We are, by 
slowly decreasing the potential, able to find the thresh- 
old at which the first domain excitation occurs and to 
demonstrate its effect on the roughness w(v). We ob- 
serve that since the domain excitation is of the order of 
the sample size, the roughness produced by that domain 
fluctuation is proportional to A. Thus there is a first- 
order jump in roughness. In contrast, a naive averaging 
of the roughness looks smooth and scales nicely. This is 
due to a scaling of the probability of a jump of the order 
A occurring at v rather than being the self-averaging be- 
havior of a typical sample. The exact numerical calcula- 
tions are supported by scaling theories based on Imry-Ma 
and large fluctuation ideas which account for the jumpy 
behavior of interfaces in a periodic potential. 

The numerical calculations are carried out using Ising 
magnets with random bonds. For a given configuration 
of bond disorder, we find the ground state interface in 
square and cubic, nearest-neighbor, spin-half, ferromag- 
netic Ising models. An interface is imposed along the 
{11} or {10} directions of a square lattice or along the 



{111} or {100} directions of a cubic lattice, by using an- 
tiperiodic boundary conditions. Periodic boundaries are 
used directions parallel to the interface unless otherwise 
mentioned. The average value of the exchange constant 
is J = 1, while the random-bond disorder is drawn from 
a uniform distribution of width S J. The periodic poten- 
tial V p — Vb[0.5 sm(2nh/X) + 0.5] is added to the random 
bond disorder, where h is to be along a direction perpen- 
dicular to the average orientation of the interface. This 
is done for the {11} and {111} cases while in the other 
orientations ({10}, {100}) we use the intrinsic lattice po- 
tential as discussed below. Notice that if A is small, the 
discrete representation of the potential will by necessity 
be rather coarse. The exact interface ground state in this 
random energy landscape is found using a mapping to 
the minimum-cut maximum-flow optimization problem 
S. We have developed a highly efficient (in both mem- 
ory and speed) implementation of the push-and-relabel 
method for the maximum flow problem 0. The exact 
ground state of a manifold in system with one million 
sites can be found in about a minute on a workstation. 



III. JUMPS BETWEEN POTENTIAL VALLEYS 

We first discuss the sensitivity of the ground state of 
the model (Q) to small variations in the amplitude Vo 
of the potential V p with wave length A. A simple scal- 
ing theory captures many aspects of this sensitivity. The 
scaling theory begins with the central limit form for the 
energy of a flat interface located at a minimum of the 
periodic potential, P\(E). If the interface is exactly flat, 
the energy fluctuations are just due to the random bond 
disorder, so that 

P I(E) * exp{-<^n, (2, 

v 7rcr I a J 

where A = L d ^ 1 is the area of the manifold and a 2 — 
2ASJ 2 is the width of Gaussian distribution. 

Consider now a system in which there are N minima in 
the periodic potential. The probability, Ljsr(E), that the 
lowest minima has energy E is, Lm{E) — NPi(E){l — 

diE)}*- 1 where, d(E) = /f^ P 1 (e) de. The difference 
in energy, g, between the lowest energy state and the next 
lowest energy state of the manifold may also be simply 
calculated. We call this difference in energy the "gap" 
and its distribution, GN{g,E) is given by, GN{g,E) = 
N(N - 1)P 1 (E)P 1 (E + g){l -d(E + g)} N ~ 2 . ' Stated 
more precisely GN(g,E) is the probability that if the 
lowest energy manifold has energy E, then the gap to 
the next lowest energy level is g. The average lowest en- 
ergy level is given by (E M ) = EL N (E) dE. This is 
not analytically tractable. However, the typical value of 
this lowest energy is estimated from aN Pi((Em)) ~ 1 
which yields, 
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J A - cr{ln(7V)} 1/2 



(3) 



To estimate the typical value of the gap, we use, 
a 2 N(N- l)Pi((E M ))Pi«.E M ) + (g)) « 1, which with 
(0) and the fact that < |(-E M )| yields, 



(.9) 



er ln(er) 

(JA-(£ M )) ~ {ln(7V)}V2 



1 H<?) 



(4) 



where cr = \/2A5J and A = L d ~ x . The gap between 
minima of the potential is thus of order l/{ln(iV)} 1//2 , 
where N ~ L^/X and L/j is the system size perpendicu- 
lar to the interface. So the separation between minima 
grows increasingly small as Lh increases. Similar extreme 
statistics problems are discussed in |jjo|| . 

Given the small gap between the metastable minima 
of the periodic potential, due to the presence of random 
bonds, we now need to find the typical change in Vo which 
can cause a level crossing in which the global ground state 
changes from one minima of the periodic potential to an- 
other |ll). The key new effect that we must control is 
the fact that the interfaces are not flat even when con- 
fined to one minimum of the periodic potential. Instead 
they have a roughness which is determined by the inter- 
play between the curvature of the periodic potential at 
its minima and the energy variations of a manifold due 
to confinement. We now develop a scaling theory for this 
phenomenon. 

First we treat the confinement effect. Consider a man- 
ifold in the presence of random bond disorder and which 
is confined in a slab of size I x L d_1 . The energy of such 
a slab is given by, 



£'!/./. = ( A 



d-l 



where L x = l 1 ^. This yields, 



c(0 



E(l,L) 



— = C\ + c 2 l~ 



where, 



x = (d-l-0)/C 



(6) 



(7) 



Notice that x is positive so that the confinement energy 
decreases as the confinement length I increases, as ex- 
pected. 

To include the effect of the confining potential, con- 
sider the behavior near a minima of the periodic potential 
to be of the form, 



V(l) = V 



A 



(8) 



where Vo = Vo /SJ, and y is a positive exponent to ensure 
that the potential is confining. For example a sinusoidal 
potential has y = 2. The behavior of a manifold in this 



confining potential and in the presence of an additive ran- 
dom bond disorder, is estimated by considering its total 
energy as a function of I [i.e. combining <Sj and (j|)], 



Ctotal = Cl 



C2 r 



Vo 



(9) 



Finding the minimum of the total energy yields the man- 
ifold roughness, 



c 2 xX« 

yV 



l/(V+a:) 



(10) 



with the energy of this optimal manifold being 



/ r v V x \ 1 /^ +;r ) 
£opt = C i+ C3 (^) , (ID 

where C3 is a constant that depends on x and y |l2| . 

Now the variation in the optimal energy with a small 
variation in Vo is given by, 



e O pt(V + 5V ) - e op i(V ) 



de pt 



5V 



(12) 



This change in energy also varies randomly from one min- 
imum of the potential to another. If the variation in the 
energy change is of order the gap found in Eq. (|]), then 
we expect the ground state location to change from one 
minimum of the potential to another. Thus we find the 
typical value of ^Vo between jumps to be found from, 



- d - lde ° pt -SVo) =(.<?}. 



dV 



(13) 



(5) Thus using Eq. (g), 



SVr 



Ojump — T d-1 



{9? ( de 



(g) 2 {x + y) f\ x v y /{x+v) 



dV J L d - 1 c 3 x V c 2 



8J % {]n{L d - 1 8J' 1 )Y f\ x V 



ln(N) 



x\r \ v/(x+v) 



SJ 



where Vo = Vo/SJ. There are several interesting features 
of this equation. Firstly note that SVojump increases loga- 
rithmically with area of the manifold On the other 
hand, the number of minima N ~ L/t/A, and SVoj ump 
decreases logarithmically with L^. The dependence of 
SVoj ump on A and on V is qualitatively as expected in 
that it increases monotonically with both of them. 

The intermittence implied by the result (|l4| ) is illus- 
trated in Fig. [l] and in Fig. ^. As a function of Vo, the 
manifold mostly stays almost unchanged in the current 
valley of minimum energy and occasionally jumps to an- 
other, new minimum of the periodic potential. A useful 
way to illustrate this intermittence as a function of Vo 
is to calculate the configurational overlap between the 
ground states as a function of V (in analogy with the 
overlap used in spin-glasses Jl3|]). The overlap, q is one 
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if two configurations are the same and zero if they have 
no bonds in common. Fig. || presents the overlap as a 
function of amplitude of the pinning potential, Vo, for 
interfaces in square and cubic lattices. The intermit- 
tent nature of periodic elastic media is clearly evident 
in these figures. Note that while the overlap and the in- 
terface roughness are intermittent, the interface energy 
(see Fig. ||) does not show any obvious signs of the jumps. 
Due to the logarithmic reduction in the gap size [Eq. (|J)], 
the interface will only sample an infinitesimal fraction 
[ln(iV)/{ln(L d_1 )} 2 ] of the available minima of the po- 
tential as we sweep v. Nevertheless a large number of 
different minima [~ ln(L^)] will be sampled by the sys- 
tem in particular if Lh is increased while the transverse 
size L is kept fixed. 



IV. ROUGHENING OF THE MANIFOLDS 

The behavior of the roughness of interfaces seen in 
Fig. H is also strongly intermittent, especially in (1 + 1)- 
dimensions. The large jumps in roughness seen in this 
figure are easily understood from the Imry-Ma argu- 
ments |l4| concerning the instability of interfaces to large 
fluctuations, as we now demonstrate for the (2 + 1)- 
dimcnsional case. 

The interface energy of a subregion a of the interface 
is, of course, also drawn from the Gaussian, P\{E) = 



eVpX 
SJ 



</tt(t 



exp 



{-^}, 



but now with standard deviation 



cr 2 = 2aSJ 2 . Some of these energy fluctuations are favor- 
able while others are unfavorable. The largest favorable 
fluctuations are found by setting ao~P\{E) rs 1, similarly 
to the extreme statistics arguments as in deriving Eq. (^), 
and it gives as the value of the energy gain, 

(E g )na{Wa)} 1/a . (15) 

A flat interface would like to take advantage of such 
large favorable energy fluctuations in adjacent minima 
of the periodic potential. However this requires having 
segments of the interface crossing the barriers in the pe- 
riodic potential. We define the barrier cost per bond to 
be e° and this is given by the integral over the barrier, 

e° = j | J^V(x)dx^ = eV - We shall use the last of 
these forms as we shall often be interested in the depen- 
dence on Vo- We consider (1 + 1) and (2 + 1) dimensional 
systems of wavelength A, length L and width B, so that 
a& — XB is the area of the part of the interface, which 
crosses the energy barrier, and a ~ LB/2 in order to 
maximize the energy gain. B = 1 is the two dimensional 
case, and B = L in the isotropic three dimensional case. 
The barrier energy cost is given by 

E b = eV a XB. (16) 

Equating Eqs. (|I1f ) and (|l6|), yields the estimate of the 
parameter values at which the first "Imry-Ma" jump in 
the manifold roughness occurs, 



'^{IniBL)}^ 



(17) 



In the (l + l)-dimensional case the logarithmic correction 
drops out, by elementary considerations. 

An "Imry-Ma" fluctuation of size a leads to a jump in 
the roughness, which is of order A x a/ A ~ A/2, since 
A = BL. We emphasize that this is the expected out- 
come in any system with fixed disorder, when Vq is var- 
ied. If B oc L, there is an exponential dependence of 
the crossover length on the parameters, for example, for 
B = L, 



L\ ~ exp 



SJ 



(18) 



an exponential dependence on v [Q. 

In Fig. |3|, we present the numerically observed behav- 
ior of the interface roughness as a function of V). We 
observe that for very large Vq the interfaces are flat and 
are confined to a minimum of the potential. For a large 
range of Vo the roughness stays the same or increases 
slowly (in three dimensions), until finally at a critical 
value a discrete jump occurs due to the "Imry-Ma" nuclc- 
ation process. This implies that the roughening process, 
as defined by the point at which the interfaces begin to 
fluctuate outside a single valley, has a first-order char- 
acter. It is seen from Fig. ^(a) that the first jump is 
oc A as expected for an extensive fluctuation. The criti- 
cal value Vo,c at which the first extensive fluctuation oc- 
curs [Fig. |(b)] follows roughly the prediction of Eq. ( |l~7j ) 
though the slope is closer to 3/4 instead of 2/3. 

The analysis of the last paragraph clearly demonstrates 
that the roughening of manifolds in periodic elastic me- 
dia is via a first order jump in roughness which is of 
order the wavelength of the periodic elastic medium. 
It is interesting to investigate whether this first order 
jump is observable in the ensemble averaged behavior. 
Scaled, ensemble-averaged plots of the manifold rough- 
ness as a function of Vq are presented in Fig. |5| for {11} 
interfaces [Fig. 0(a)] and for {111} oriented interfaces 
[Fig. ||(b)]. These plots scale quite nicely with the char- 
acteristic length and roughness suggested by Eqs. ( |l7| ) 
and (|isj). In the two dimensional case, there is also a 
clear indication of the first order character of the transi- 
tion. The three dimensional data gives little indication 
of the first order jump in roughness and underscores the 
problems with a naive averaging of the data. However, 
we do not have any clear explanation, why the roughness 
values in the plateau before the jump can be collapsed 
with the same prefactor as in the asymptotic roughness 
in the {111} case, but not in the {11} case. 
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V. PERIODICITY DUE TO THE LATTICE 



It is of interest to see if the first order character of the 
roughening of PEM extends to manifolds in the {10} and 
{100} directions. In these directions, the lattice itself in- 
troduces a periodicity, which for example is the origin 
of the thermal roughening transition in lattice models 
in three dimensions. Thus we do not need to introduce 
an extra periodic potential, and instead we just study 
the roughness of these manifolds as a function of disor- 
der. We have studied the roughness of {100} manifolds 
as a function of disorder before and in those studies we 

In the light of the un- 



ensemble averaged the data 15 



derstanding develop above, we have revisited this prob- 
lem and find that the typical behavior in both the {10} 
and {100} problems is very similar to that suggested by 
the PEM model. That is, in a large sample the system 
roughens via a first order jump in the roughness due to 
an extensive fluctuation. The behavior of one sample 
as a function of disorder is presented in Fig. §(a). The 
probability distributions of the roughness for several L 
are presented in Fig. 0(b) in which we observe how one 
can pass through a co- existence region with both flat and 
rough samples as L is varied. The intermittent behavior 
typical of PEM is evident in Fig. 0(a), but is obscured 
by the averaging in Fig. 0(a). Though the jump transi- 
tion from the flat phase to the algebraically rough phase 
occurs in both the periodic elastic model in the {111} 
direction and for interfaces in the {100} direction, there 
is an important difference in the behavior of these mod- 
els [compare Figs. |(b) and 0(a)]. In the PEM model 
in the {111} direction, there is a pronounced plateau in 
the roughness due to the saturation of wandering within 
one well [Fig. |](b)]. In contrast, in the {100} direction, 
the interface remains flat until the transition to the alge- 
braically rough phase [see Fig. 0(a)]. The extent of the 
plateau region can be tuned in the PEM model by vary- 
ing the shape of the potential near the minimum and by 
varying the wavelength. We have also carried out calcula- 
tions for the case of dilution disorder [Fig. 0(b)] and find 
similar behavior, with the averaged behavior presented in 
Fig. 0(a). With dilution disorder the pronounced plateau 
is not due to any roughening inside a valley, but because 
of rare "bumps" , whose occurrence is due to the Poisso- 
nian statistics of diluted bonds. The averaged data scales 
quite well with (SJ/J) 2 = p(l - p) J 2 /{pJ) 2 = (l-p)/p, 
where J = 1 and the variance of the binomial distribu- 
tion var = std 2 = p(l — p)J 2 with the corresponding 
mean p J, and p is the occupation probability of a bond. 
Thus we find, in contrast to our earlier conclusions from 
similar data, that at large enough length scales interfaces 
in the {100} orientation are algebraically rough and are 
consistent with the PEM model. 

A further important feature of the large fluctua- 
tion character of the roughening transition is that it is 
strongly dependent on the boundary conditions. This 



is illustrated in Figs. 0(b) and 0(a) in which the rough- 
ness is depicted as a function of the amplitude of the 
disorder for both periodic and free boundaries and with 
the same disorder configuration. The threshold value of 
Vb at which the first order jump in roughening occurs 
is typically smaller for the case of periodic boundaries. 
Large fluctuations can take advantage of the boundary 
to reduce the cost of crossing the energy barrier. This 
sensitivity to boundary conditions is a hallmark of the 
large fluctuation effects discussed here. 



VI. CONCLUSIONS 

To conclude, we have discussed the roughening of elas- 
tic manifolds in the presence of a competition between 
bulk randomness and a confining periodic potential. We 
have concentrated on the two- and three-dimensional 
cases which are well-known to have, asymptotically, an 
algebraic roughness scaling. A study of the system-by- 
system behavior reveals however a much richer scenario 
in which each manifold makes intermittent jumps, finally 
culminating in a first-order change in its roughness. This 
process is also important since it is related to the asymp- 
totic scaling of the roughness. Recent experiments on the 
creep of (l+l)-dimensional systems [Q have shown that 
scaling arguments of activation energy barriers can match 
real systems, using predictions based on rough manifolds. 
The time scales also depend crucially on the actual am- 
plitude which is set in our picture by the roughening 
transition. 

Also, the intcrmittcncc in the early stages would merit 
experimental consideration. Such jumps in the mean lo- 
cation of the interface, could be studied in the asymp- 
totic rough regime. In an independent study we have 
pointed out this mechanism for both fracture surfaces 
arising from random fuse networks, and from yield sur- 
faces of perfectly plastic media which are equivalent to 
the minimum energy surfaces studied here jlq ]. 

The focus of renormalization group and variational 
calculations in this problem has been dimensions d = 
(D + 1) > 4, since there one encounters two asymptotic 
regimes separated by a transition. Of the two phenom- 
ena discussed here at least the intermittent jumps in the 
center of mass location of the interface should persist in 
that case. 
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FIG. 1. Interface configurations in (1 + l)-dimensions for 
various Vb- In this calculation the disorder configuration and 
wavelength (A = 16) are fixed at S J — 1. As Vb is varied, the 
interface jumps between the minima of the periodic potential. 
The solid lines denote the position of the largest values of the 
sinusoidal periodic potential V p . The lattice size is 160 x 160 
and the interfaces are oriented along the {11} direction. Note 
that the disorder is exactly the same for each value of Vb- 
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FIG. 2. The overlap q = L~ (d_1) ^(/i, 1 - h'f) between 
ground states as the amplitude of the potential Vo is var- 
ied (SJ = 1). As Vo is decreased, we calculate the overlap 
between the interface configuration at one value of Vo (de- 
scribed by {hi}) and the interface at the next value of Vo 
(described by {h 2 }). The corresponding mean heights (h) 
are shown in the insets. The calculations were carried out 
as for Fig. |lj, however we used 300 different values of Vo with 
AVo = 10~ for the same realisation of disorder and the wave- 
length A = 4. (a) Two dimensional case, with the system 
size L x Lh — 1024 x 1025. (b) (2 + l)-dimensional inter- 
faces oriented along the {111} direction for lattices of size 
L 2 x L h = 100 2 x 129. 
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FIG. 3. The interface width [w 2 = L~ (d - 1) J2i( h * ~ ( h )) 2 ] 
and the total energy as a function of Vo for A = 4 and SJ = 1. 
The results are for a fixed disorder configuration and from the 
same calculations as Fig. ^| (a) (1 + l)-dimensional system, 
(b) (2 + l)-dimensional system. The systems with free and 
periodic boundaries have the same realization of randomness. 
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FIG. 4. (a) Average size of the first jump in roughness Aw, 
when Vq — Vo,o normalized using A and calculated as the dif- 
ference between roughness values just after a jump and before 
that, as a function of the volume of the systems. We have 
carried out simulations for a strips of dimension L — 1000, 
B = 1 — 64, and Lh ~ 5A, for various values of A. The num- 
ber of realizations is 100. (b) Average value of the amplitude 
of the potential Vo = Vo, c at which the large-scale "Imry-Ma" 
fluctuation occurs (5 J =1). The data is from the same simu- 
lations as in (a) for (2+ 1) dimensional case, i.e. B > 1. The 
results are scaled using the prediction ( |l7| ) . 
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FIG. 5. Scaled roughness of interfaces oriented in {11} and 
{111} directions, for various values of Vo and L. (a) {11} 
oriented systems with A = 16, SJ = 1, and system sizes 



L = 20 A 



1280 . The number of realizations is 200 for 



each system size and Vo- The solid line corresponds the slope 
£ = 2/3. (b) {111} oriented systems with A = 4 and system 
sizes L 3 = 10 3 - 90 3 . The number of realizations is 200 for 
each system size, SJ and Vo- The solid line corresponds the 
slope £ = 0.42, while the dotted line is £ = 0.36. 
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FIG. 6. Behavior of the roughness of interfaces oriented 
in the {100} direction, (a) The intermittence of a single re- 
alization as a function of the amplitude of uniform disorder 
SJ. The disorder configuration is the same (with both free 
and periodic boundaries), but the ratio 6 J/ J is slowly in- 
creased in steps of 0.01. The system size is L 3 = 100 3 . (b) 
The histograms of the roughness values w for system sizes 
L 2 x L h = 50 3 . . . 200 2 x 100. The peak of the distributions 
jumps from w ~ to w ~ 0.5, when the system size increases. 
The number of realizations is 500 for smaller system sizes and 
200 for L 2 x L h = 200 2 x 100. SJ/J = 0.9. 
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FIG. 7. (a) Scaled roughness for continuum disorder. The 
system sizes ranges from L 2 x Lh = 4 3 to 200 2 x 100. 
The number of realizations ranges from 500 for system sizes 
L 2 x L h < 140 2 x 100 to 200 for the larger ones, (b) Scaled 
roughness for dilution type of disorder. The system sizes 
ranges from L 2 x Lh — 4 3 to 200 2 x 100 (and even up to 
400 3 for p — 0.90). The number of realizations ranges from 
500 for system sizes L 2 x L h < 140 2 x 100 to 200 for the larger 
ones (with the exception of larger system sizes for p — 0.90). 
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